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We introduce a one-dimensional model of Bose-Einstein condensates (BECs), com- 
bining the double-well potential, which is a well-known setting for the onset of 
spontaneous-symmetry-breaking (SSB) effects, and time-periodic modulation of the 
nonlinearity, which may be implemented by means of the Feshbach-resonance- 
management (FRM) technique. Both cases of the nonlinearity which is repulsive 
or attractive on the average are considered. In the former case, the main effect pro- 
duced by the application of the FRM is spontaneous self-trapping of the condensate 
in either of the two potential wells in parameter regimes where it would remain un- 
trapped in the absence of the management. In the weakly nonlinear regime, the 
frequency of intrinsic oscillations in the FRM-induced trapped state is very close to 
half the FRM frequency, suggesting that the effect is accounted for by a paramet- 
ric resonance. In the case of the attractive nonlinearity, the FRM-induced effect is 
the opposite, i.e., enforced detrapping of a state which is self-trapped in its unman- 
aged form. In the latter case, the frequency of oscillations of the untrapped mode is 
close to a quarter of the driving frequency, suggesting that a higher-order parametric 
resonance may account for this effect. 

PACS numbers: 03.75.Lm: 42.65.Wi; 03.75.Lm 
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Symmetric double-well potentials (DWPs) play an important role in many 
physical situations dominated by self-focusing and defocusing nonlinearities. It 
is well known that, in quantum mechanics, the DWP gives rise to alternating 
symmetric and antisymmetric bound states. However, the nonlinearity breaks 
the symmetry, giving rise to asymmetric bound states. This effect of the spon- 
taneous symmetry breaking (SSB) has drawn a great deal of attention in the- 
oretical and experimental studies of nonlinear physical media, especially in the 
context of Bose-Einstein condensation (EEC). In this work, we aim to extend 
the analysis of the SSB by introducing a BEC model which combines the DWP 
and the time-periodic modulation of the strength of the nonlinearity. The latter 
is another important tool used for the control of the nonlinear dynamics in BEC, 
by means of the physical effect known as the Feshbach resonance (controlled by 
a time-periodic external field, in that case). We consider the interplay of the 
DWP and time- modulated Feshbach resonance {Feshbach-resonance management, 
FRM) for both signs of the nonlinearity, repulsive and attractive. In these two 
cases, new dynamical effects are predicted by means of a systematic numerical 
analysis. For the repulsive (on the average) nonlinearity, the application of the 
FRM induces self-trapping of the condensate in one of the two symmetric poten- 
tial wells, in the case when the free condensate would be untrapped. In such 
an FRM-induced trapped state, for low density cases, the frequency of intrinsic 
oscillations is found to be very close to half the underlying FRM driving fre- 
quency, which suggests that the trapping is induced by a parametric resonance. 
In the opposite case of the attractive (on the average) nonlinearity, the effect of 
the FRM is opposite too, namely it is detrapping of states which were self-trapped 
in the absence of the FRM, with the frequency of oscillations of the detrapped 
condensate close to a quarter of the driving frequency for low density cases too. 
A higher-order parametric resonance may plausibly account for the latter effect. 

I. INTRODUCTION 

The past two decades have witnessed remarkable progress in experimental and theoretical 
studies of Bose-Einstein condensates (BECs)^'^. Nonlinear matter-wave structures have 
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been one of the significant subjects considered in this context^"-. From the theoretical 
viewpoint, the emergence of such structures can be understood in the framework of the well- 
estabhshed mean-field approximation, based on the Gross-Pitaevskii equation. In particular, 
for attractive (repulsive) interatomic interactions, this equation predicts the existence of 
bright (dark) matter- wave solitons, which have been observed in a series of experiments, see 
Refs." and^'"— , respectively. Bright gap solitons are also possible in BEG with repulsive 
interactions loaded into an optical lattice^^. 

An important element in the studies of BEC is the variety of external potentials that may 
be used to confine the ultracold atomic gases. The most typical forms of such trapping po- 
tentials are the harmonic-oscillator traps and optical lattices, i.e., periodic potentials created 
by the interference of counter-propagating laser beams. A similar situation is relevant also 
in the context of nonlinear optics, where the basic model, namely the nonlinear Schrodinger 
equation, may incorporate harmonic or periodic potentials describing graded-index waveg- 
uides and periodic waveguiding arrays, respectively^. A combination of harmonic and 
periodic traps may be used to create a double- well potential (DWP), which has drawn a 
great deal of attention. BECs confined in DWPs were studied experimentally in Refs.— i^, 
and many theoretical works addressed this setting^^i"— . A fundamental effect studied, in di- 
verse forms, in these works, is the spontaneous symmetry breaking (SSB) of the population 
of atoms in the two symmetric wells. A related work was done in nonlinear optics, where 
twin-core self-guided laser beams in Kerr mediai^, optically induced dual-core waveguiding 
structures in photorefractive crystals2£, and trapped light beams in an annular core of an 
optical fiber— among others, also lead to manifestations of phenomenology associated with 
DWPs. 

In addition to the studies of the SSB in the static DWP settings, the analysis was extended 
to diverse situations with the DWP subject to "management", i.e., its parameters, such the 
depth of the well or height of the barrier between them, were made periodically varying 
functions of time'^^. In fact, in all those works the analysis was based on a two- mode 
approximation, which replaces the underlying Gross-Pitaevskii equation (GPE) by a system 
of two coupled ordinary differential equations (DDEs). 

Apart from the variety in the shapes of the potentials, another important tool for manip- 
ulations of BEG is provided by magnetically^*^ or optically^! induced Feshbach resonance, 
which makes it possible to control the effective nonlinearity in the condensate. The lat- 
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ter possibility has given rise to many theoretical and experimental studies. A well-known 
example is the formation of bright matter-wave solitons and soliton trains in ^Li^i^ and 
^^Rb^ condensates, by switching the interatomic interactions from repulsive to attractive. 
Many theoretical works studied the BEC dynamics under temporal and/or spatial mod- 
ulation of the nonlinearity. In particular, the application of such a "Feshbach resonance 
management" (FRM) technique in the temporal domain can be used to stabilize attractive 
higher-dimensional BEC against collapse^*^, and to create robust matter-wave breathers 
in the effectively one- dimensional (ID) condensate^. On the other hand, the so-called "col- 
lisionally inhomogeneous condensates", controlled by the spatially modulated nonlinearity, 
have been predicted to support a variety of new effects, such as the adiabatic compression of 
matter-waves^!^, Bloch oscillations of matter-wave solitons^, atomic soliton emission and 
atom lasers^, dynamical trapping of matter-wave solitons^i^, enhancement of transmis- 
sivity of matter-waves through barriers^*^, stable condensates exhibiting both attractive 
and repulsive interatomic interactions^"^, the delocalization transition^, SSB in a non- 
linear double-well pseudopotential^, the competition between incommensurable linear and 
nonlinear lattices^, and others. 

In this work we combine the two above-mentioned settings, namely the DWP and a time- 
modulated nonlinearity, in the ID geometry, with the objective to control the SSB in the 
double-well setting by means of the FRM. Our model is based on the corresponding GPE, 
see Eqs. ([1]) and ([2]) below. We consider the nonlinearity which may either be repulsive 
or attractive on the average. In the case of the repulsive interactions, we find that the use 
of the FRM results in the onset of spontaneous self-trapping of the BEC in one of the two 
potential wells, while without the application of the FRM atoms oscillate between the wells, 
i.e., the condensate is untrapped. We also find that, when the self-trapping occurs, the 
trapped BEC undergoes small- amplitude intrinsic oscillations at a frequency equal to half 
the nonlinearity-modulation frequency, which indicates the role of a parametric resonance in 
this case, when the densities are sufficiently low. In the case of the attractive nonlinearity, 
the FRM-induced effect turns out to be the opposite: while the unmanaged condensate is 
spontaneously trapped in one of the two potential wells, the application of the FRM leads 
to its detrapping (dynamical symmetry restoration), with the density oscillations between 
the wells at a frequency which is almost exactly a quarter of the driving frequency for low 
density cases. Thus, for either nonlinearity (repulsive or attractive), our results show that 
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the FRM can control the BEC dynamics in DWPs, by inducing trapping/detrapping and 
intrinsic oscillations, which do not take place in the absence of the FRM. Our numerical 
results, obtained in the framework of the GPE model, are also corroborated by a semi- 
analytical approximation. In particular, we extend the analytical approximation presented 
in Refs.— to the FRM-driven setting under consideration: we adopt a Galerkin-type 
expansion to describe the evolution of the wavefunctions at each well of the DWP. This 
leads to a non-autonomous system of ODEs, which are solved numerically. In all cases, the 
ODE result is compared to the one obtained by the GPE model and, in most cases, they 
are found to be in fairly good agreement. 



It is relevant to mention that a similar physical model was recently introduced in Ref.— , 
where it was postulated that, under the action of the FRM driving, the DWP was described 
by a system of linearly coupled ODEs corresponding to the above-mentioned two-mode 
approximation. However, dynamical regimes studied in Ref.™ were fairly different from 
those considered here. The analysis reported in that work demonstrated an FRM-induced 
transition from untrapped oscillations to a trapped state, but this was observed in the high- 
frequency regime, and was explained by means of an averaging method. In the present work, 
the transition to the trapping is clearly accounted for by a parametric resonance at moderate 
values of the FRM frequency u. of the underlying GPE. On the other hand, the ODE 
system used in Ref.— cannot explain the parametric resonance (clearly, a more sophisticated 
approximation is needed to capture the resonant mechanism). Another difference is that 
we observe the FRM-induced trapping and detrapping effects (the latter was not reported 
in Ref.—) in the regime of relatively weak FRM, when the time- dependent nonlinearity 
coefficient, g{t), does not change its sign. In Ref.—, the opposite regime of the strong 
management was considered. 



The rest of the paper is structured as follows. In Section II we introduce the model, 
and present the results for the repulsive and attractive nonlinearity in Sections III and IV, 
respectively. Section V concludes the paper. 
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II. THE MODEL AND ITS CONSIDERATION 



A. The model 

In the mean-field approximation, the combination of the DWP and time-varying non- 
linearity is described by the GPE for the macroscopic wavefunction ip {x,t), written in the 
scaled form: 

^^l^t = - (1/2) + U{x)i; + gitM^^ - fi^. (1) 
Here is the chemical potential, while the DWP is taken in the form of^^i*^ 

U{x) = (1/2) fi^a;^ + l/osech^^ (x/w) , (2) 

where Q is the normalized strength of the harmonic trap, while Vq and w denote, respectively, 
the amplitude and width of the barrier used (in combination with the parabolic trap) to 
create the DWP. Finally, the FRM modulation function g{t) is taken in the customary 
form^, 

9{t) = go + 9i sin (ut) . (3) 

The same model may be interpreted in terms of nonlinear optics, with t being the prop- 
agation distance, /x the propagation constant, and g(t) representing a periodic modulation 
of the Kerr coefficient^. In that case, potential ([2]) defines two waveguiding channels in a 
planar medium. 

We aim to study localized FRM-driven states in the ID setting described by Eqs. ([1]) 
and ([2]). In particular, we numerically solve the above model with initial conditions of the 
form of either an antisymmetric (with respect to the two wells) initial state and go > (the 
self-repulsive nonlinearity, on the average), or a symmetric state and gQ < (i.e., the non- 
linearity which is self-attractive, on the average). This choice stems from the fact that both 
configurations are well known to feature the SSB, in the absence of the "management"—"—. 

Numerical results, which are reported below, were obtained by means of the split-step 
Fourier method. We use the following values of parameters which adequately represent the 
generic situation: the effective parabolic trap strength Q assumes value 0.1, the chemical 
potential, /i, ranges from 0.05 to 0.5, the barrier's width is w = 0.5, and its height is Vq = 1. 
As concerns the nonlinearity-modulation function, g{t), the "dc value" go is taken as ±1 
for the repulsive or attractive interatomic interactions, respectively, while the modulation 
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strength, gi, varies in the interval of (0,1.5). Finally, the driving frequency takes values 
in the interval of < u; < 1. In the context of BEC, such parameters may correspond 
to the 23Na or a ''Li condensate, for the repulsive {qq = +1) and attractive {qq = — 1) 
nonlinearity, respectively, confined in the trap with longitudinal and transverse confining 
frequencies 2tt x 10 Hz and 27r x 100 Hz. In these two cases, the numbers of ^^Na and ^Li 
atoms are, respectively, ~ 500 and 1000. The corresponding peak values of the normalized 
peak densities (at the DWP minima) are (I^Hjuax ~ ^-^ ^■^'^ repulsive case, or 

0.13 for the case of the attraction, for the respective chemical potential is 0.22, 0.3 or 0.05. 

B. Semi- analytical approach 

The spectrum of the Schrodinger equation Eq. ([1]) for g(t) =0, consists of a ground state, 
ipo{x)j and excited states, ipilx) {I > 1). In the weakly nonlinear regime of the full problem, 
using a Galerkin-type approach, we expand ip{x,t) as^, 

^{x, t) = Co{t)Mx) + Ci(t)V^i(x) + ■ • ■ , (4) 

and truncate the expansion, keeping solely the first two modes; here Co,i(t) are unknown time- 
dependent complex prefactors. Once again, it is worth noticing that such an approximation 
(involving the truncation of higher-order modes and the spatio-temporal factorization of 
the wavefunction) , is expected to be quite useful for a weakly nonlinear analysis, i.e., for a 
sufficiently small norm (or, physically speaking, number of atoms) of the solution. 

Substituting Eq. (jlj) into Eq. ([T]), and projecting the result onto the corresponding eigen- 
modes, we obtain the following ordinary differential equations (ODEs), for the coefficients 
of the projection onto ipo and ipi, respectively: 

ico = {uo - /i)co - (5-0 + fi'i(t))Ao|copco - (5-0 + fi'i(t))5(2|ci|^co + c^cq), (5) 
ici = {uji - /i)ci - {go + gi{t))Ai\ci\'^ci - (go + gi{t))B{2\co\'^ci + CgCi). (6) 

In Eqs. ([5])- ([6]), dots denote time derivatives, the overbar stands for the complex conjugate, 
Wo,! are eigenvalues corresponding to the eigenstates ipo,!, while Aq = f ip^dx, Ai = f ipfdx, 
and B = J tp^tpfclx are constants. Additional overlapping integrals, viz., J ipQipfclx and 
J ipiip^dx, which formally appear in the course of the derivation of the ODE system, vanish 
due to the opposite parities of the real wavefunctions ipQ and ipi, in the framework of the 
underlying linear Schrodinger problem with the symmetric potential. 
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We now use amplitude-phase variables, Cj = pje^'^'^ , j = 0,1 (the amplitudes pj and 
phases (pj are assumed to be real), to derive from the ODEs ([5])-([6]) a set of four equations. 
Introducing function ip = (pi — (po, we obtain the equations for po and (pQ-. 

Po = {go + gi{t))Bp\pQ sin(2¥?), (7) 
00 = (/X - uo) - (^0 + 9iit))Aopl - 2{go + gi{t))Bpl - {go + gi{t))Bpl cos(2</;), (8) 

while the equations for pi, (pi are 

Pi = -{% + gi{t))Bplpi sin(2(^), (9) 
^i = {li- ui) - {go + gi{t))Aopl - 2{go + gi{t))Bpl - {go + gi{t))Bpl cos(2(/^), (10) 

The conservation of the norm N = J {ipl'^dx implies Po + Pi = N. Finally, subtracting 
Eqs. ([8]) and ffTOl) . we obtain: 

^ = -Au + {go + gi{t)){Aopl - Aipl) - {go + gi{t))B{2 + cos{2^)){pl - p?). (11) 

where Au = ui — Uo- Equations ([7]), Q) and (ITTl) constitute a non-autonomous dynamical 
system, which we solve numerically. This way, from (po, (po) and (pi, (pi) we can respectively 
find Co{t) and Ci{t) and, thus, obtain the wavefunction profile as per the approximation of 
Eq. (jll). For a given set of parameter values, the outcome of this calculation will be compared 
to the outcome of the direct numerical integration of Eq. ([T]). 

Below we consider, at first, the case of the repulsive nonlinearity, and separately an- 
alyze the situations corresponding to "small" and "large" initial peak densities, with the 
corresponding values (I'^Hmax = 0-1 0.17. 

III. THE REPULSIVE NONLINEARITY 

A. "Small" initial values of the peak density 

In the case of the self-repulsive nonlinearity, we start by considering the case of the 
"unmanaged" condensate {gi = 0), with parameters taken as in Ref.— , viz., go = 1, = 0.1, 
Vo = I, w = 0.5, and initial value (iV^Hmax ~ Moreover, the parameters involved in Eqs. 
(CD, (E]) and (HH) are found to be Ao = 0.09078, Ai = 0.09502, B = 0.08964, coo = 0.13282 
and Ui = 0.15571. In the left panel of Fig. ([1]) we show a spatiotemporal contour plot 
of the evolution of the density, while the right panel of the same figure displays the time 
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evolution of the relative difference in the atomic population between the left and the right 
wells, defined as 



n{t) 



(12) 



Note that in the right panel of Fig. ([1]) we show n{t) as obtained from the GPE (solid line) 
and the ODEs, Eqs. ([7]), ([9]) and f|TT]) (dashed line); the agreement between the two is very 
good. A variable similar to n{t) was also used in the recent work^^, which analyzed the 
present setting in the framework of the ODEs of a two-mode reduction (see e.g. Ref.— ). 
The free oscillations of n{t) in the untrapped state feature a well-defined frequency, which 
can be estimated from the right panel of Fig. [T]as cJqsc ~ 0.021. Moreover, from the value 
of n{t) one can straightforwardly estimate the fraction of the total number of atoms which 
oscillates between the two wells and which is greater as the value of n{t) approaches the ±1. 
In the case of Fig. [T]the maximum value of \n{t)\ is 0.84. 
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FIG. 1. (Color online) Left panel: The spatiotemporal evolution of the density in the double-well 
potential. Right panel: the evolution of the relative difference between the numbers of atoms in 
the left and right wells, n{t). The solid line represents the results from the integration of the 
GPE, while the dashed one the results obtained from the ODEs. In this case, the nonlinearity is 
repulsive, go = 1, with gi = and initial value (|V'l^)j„ax ~ ^'^^ while parameters of the potential 
are = 0.1, Vq = 1 and w = 0.5. 



The next step is to vary the FRM frequency, u, at different fixed values of the FRM 
strength, gi. Starting with gi = 0.2, we varied u from to 1 with step 1 x 10~^. The result 
was that, in the intervals of < u < 0.0793 and u > 0.1677, the oscillations of n{t) remain 
periodic, like in Fig. [H with almost the same frequency as at gi = 0, i.e. Uosc ~ 0.021, 
and with the zero average value of n{t), i.e., the managed system remains in the untrapped 



state. 

However, the example displayed in Fig. |2]for uj = 0.11 demonstrates that the situation 
is completely different for 

0.0793 <uj< 0.1677, i.e., 4u;osc < oo < Setose , (13) 

when the the system gets trapped in either of the two wells, thus exhibiting an FRM-induced 
SSB. Actually, for all the values of u from interval (fT3|) . n{t) does not change its sign, while 
the largest values of \n{t)\ are above 0.5, hence the observed trapping is quite strong. In the 
example shown in Fig. [21 the observed frequency of the trapped oscillations is very close 
to half the FRM frequency: Wtrap = 0.054 ^ uj/2 (note that Utrap is very different from 
the above-mentioned frequency of the oscillations in the untrapped state, Uosc ~ 0.021). 
Examining the trapped states for other values of u from interval (fT3!) . we have concluded 
that relation Utrap ~ w/2 remains valid as long as the FRM-induced trapping holds, as seen 
in the left panel of Fig. [3l which displays the dependence between the values of u and Uosc 
in the domain of 4a;osc < uj < Suosc- This relation clearly suggests that the SSB trapping 
induced by the FRM is a result of a parametric resonance— . To the best of our knowledge, 
this is the first example of the resonant SSB in the DWP model, induced by the time-periodic 
drive. It is relevant to mention that the parametric resonance also plays a dominant role 
in effects induced by the management in the form of the time-periodic modulation of the 
strength of the single-well trapping potential^. 




t t 

FIG. 2. (Color online) The same as in Fig. [H but in the presence of FRM with gi = 0.2, lo = 0.11, 
and initial vahie (IV'P)juax ~ right panel, the solid and dashed lines represent the results 

obtained from the integration of the GPE and ODEs, respectively. Other parameters are as in Fig. 

m 
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Next, increasing the value of the FRM amphtude to gi = 0.4, the results are as follows: 
for < < 0.0863 and u > 0.1581, the condensate remains untrapped, with n{t) featuring 
the oscillation frequency which is virtually identical to that observed in the absence of the 
FRM (i.e. Uosc ~ 0.021). On the other hand, in the interval of 0.0863 < to < 0.1581, which 
is almost identical to 4wosc < w < 8u;osc, cf. its approximate counterpart in Eq. (IT^ . the 
oscillations of n(t) demonstrate the SSB leading to a state trapped in one potential well 
with maximum values of \n{t)\ greater than 0.5. In this case too, we have examined the 
relation between the FRM frequency, u, and the frequency of the trapped oscillations, entrap- 
The obtained results corroborate that the above-mentioned parametric-resonance relation, 
f^trap ~ Ci;/2, remains valid, as seen in the right panel of Fig. |3J 




FIG. 3. (Color online) The relation between the FRM frequency, u, and frequency Wtrap of the 
oscillations of the condensate spontaneously trapped in one well under the action of the FRM, 
in the case of the repulsive nonlinearity. Solid, and dotted lines correspond, respectively, to the 
result obtained from the GPE, and the estimation of Wtrap = Left panel: gi = 0.2 and 

0. 08 < oj < 0.1677. Right panel: gi = 0.4 and 0.09 < uj < 0.15. Other parameters are as in Figs. 
Elandll 

Increasing the value of gi, we observed qualitatively the same results up to gi = 0.9, 

1. e., as long as the FRM amplitude, gi, remained smaller than the average value of the 
repulsive-nonlinearity coefficient, go = 1. For gi > go, when the total nonlinearity coefficient 
periodically changes its sign, the situation is qualitatively different. In particular, for gi = 
1.2, the SSB (trapping in one well) occurs at 0.002 < u < 0.011 and 0.088 < a; < 0.097 
[see an example in the left panels of Fig. HI for w = 0.005]. At 0.011 < u < 0.088, 
the condensate remains untrapped [see middle panels of Fig. HJ for u = 0.025], while at 
u > 0.097 the oscillations of n{t) are irregular, featuring neither symmetric oscillations 
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nor the self-trapping [see right panels of Fig. HJ for uj = 0.15]. The parametric- resonance 
relation, cutrap ~ i^/2, is vahd only in the interval of 0.088 < a; < 0.097. 

Similar results were obtained for gi = 1.5. In this case, the SSB trapping occurs in the 
interval of 0.001 < u < 0.012 [see left panels of Fig. [5l for w = 0.005]. At u > 0.012, 
n{t) again features the irregular evolution, see an example in the right panels of Fig. [Sj 
for u = 0.055. Note that in these cases, corresponding to relatively large value of gi, the 
Galerkin-type approximation for n{t) fails to follow the result obtained in the framework of 
the GPE. This may be natural to expect as the large strength of the nonlinearity renders 
relevant the inclusion of additional modes in the description of the BEG dynamics. 
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FIG. 4. (Color online) The top panels show the density contour plots and the bottom panels the 
oscillations of the population imbalance between the two potential wells, n{t), in the case of the 
repulsive nonlinearity with go = 1, gi = 1.2, and initial value of (I'i/'P) = 0.1. Left panel: 

^ ' 1X1. 

u} = 0.005; middle panel: uj = 0.025; right panel: oj = 0.15. The parameters of the potential are 
the same as in the previous figures. In the bottom panels, the solid and dashed lines represent the 
results obtained from the integration of the GPE and ODEs, respectively. 

B. "Large" initial values of the peak density 

Here, we extend the analysis, increasing the initial values to (|'i/'P)jnax ~ 0.17, while all 
the other DWP parameters are the same as above (i.e. = 0.1 and Vq = 1). For gi = 
(no management), the results are shown in Fig. |6l In this case, the frequency of the free 
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FIG. 5. (Color online) The evolution of the density contour plots (top) and of n(i) (bottom) in the 
case of the repulsive nonlinearity with go = 1, gi = 1.5, and initial value (iV'l^)n^ax ~ ^^^^ and 
right panels pertain to a; = 0.005 and oj = 0.055, respectively. The potential is the same as in the 
cases considered above. Additionally, In the bottom panels, the solid and dashed lines represent 
the results obtained from the integration of the GPE and ODEs, respectively. 

oscillations of the condensate between the two wells is Wqsc ~ 0.031. The maximum value 
of the relative population imbalance, n{t), is 0.602, and for larger values of the initial peak 
density this value decreases. 

Next, we study the behavior of n{t) for different values of gi and u. For gi = 0.2, the result 
is that the condensate remains untrapped at < < 0.06341 and u > 0.1009 [similar to 
the situation displayed in Fig. |6], featuring the same oscillation frequency as in the absence 
of the management, Wqsc ~ 0.031. For 0.06341 < u < 0.1009 (i.e., 2u;osc ^ ^ Swqsc), the 
evolution of n{t) suggests trapping in one well, demonstrating the SSB effect. This result is 
shown in Fig. [7]for = 0.08. Notice that, as observed in the right panels of Figs. [6] and [71 
the agreement between the ODEs and the GPE results is fairly good. 

Here, we should recall that, in the previous case, with the smaller density peaks, the 

13 




I 

-I 
-I 

I 



200 400 600 800 1000 
t 



0.25 


Ir 


0.2 


0.5 


0.15 


S 0- 


0.1 




0.05 


-0.5 

-i 




200 



400 600 
t 



800 1000 



FIG. 6. (Color online) The same as in Fig. [H but for the initial value (|'0P)j„ax ~ 0.17. In the 
right panel, the solid and dashed lines represent the results obtained from the integration of the 
GPE and ODEs, respectively. The parameters of the potential are the same as in the previous 
case. 



observed frequency of the oscillations of the trapped BEC was following the parametric- 
resonance relation, Wtrap = uj/2. As the initial density peak is getting larger, and, in partic- 
ular, in the present case it corresponds to (I'^AHmax ~ 0.17, this rule is no longer valid. In 
particular, the frequency of the trapped oscillations in the case of u = 0.08 is 0.0604. 
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FIG. 7. (Color online) The same as in Fig. [21 but for gi = 0.2, oj = 0.08, and initial value 
(l^l^)max ~ 0.17. In the right panel, the solid and dashed lines represent the results obtained from 
the integration of the GPE and ODEs, respectively. The other parameters are as in Fig. [2j 



Next, we increase the FRM amplitude, gi. For gi = 0.4, we have found that, at < 
u < 0.0598 and u > 0.0904, the trajectory of n(t) remains untrapped, with the same 
oscillation frequency as in the absence of the management. At 0.0598 < u < 0.0904 (i.e., 
ScUosc ^ ^ ^ Scjosc), the oscillations of n{t) feature the self-trapping in one well. 
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Qualitatively the same behavior is observed for gi < 0.95. At greater values of gi, we did 
not observe any self-trapping. More specifically, we examined the same two values of the 
FRM amplitude as above, i.e., gi = 1.2 and gi = 1.5. For gi = 1.2, we observed untrapped 
periodic oscillations for < a; < 0.0561 [see the left panel of Fig. [H]for u; = 0.01], while for 
u > 0.02 the population imbalance, n{t), features an irregular evolution, see the right panel 
of Fig. [8]for u = 0.1. Notice that, again, the agreement between the ODEs and GPE result 
is less good for this case of large nonlinearity due to the apparent involvement of higher 
modes in the dynamics. 

Similar results were obtained for gi = 1.5. In this case, untrapped periodic oscillations 
occur at < < 0.001, resembling the situation displayed in the left panel of Fig. IH while, 
at w > 0.001, irregular evolution of n{t) is again observed, similar to that in the right panel 
of Fig. El 
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t t 

FIG. 8. (Color online) Trajectories of n{t) in the case of the repulsive nonlinearity with gQ = 1, 
gi = 1.2, and initial value (|''/'l^)n^ax ~ 0-17. Left and right panels pertain, respectively, to w = 0.01 
and oj = 0.1 with the same potential parameters as above. Moreover, the solid and dashed lines 
represent the results obtained from the integration of the GPE and ODEs, respectively. 

In Table I, we summarize intervals of gi/u in which the trapped oscillations take place, 
in the cases of smaller and larger initial values of the peak density {\ip\'^)max [i-e., 0.1 and 
0.17], for gi = 0.2 and 0.4. In these cases, we observe that the trapped oscillations appear 
at values of gi/u) similar to those presented in Ref.— , in the framework of the two- mode 
approximation. However, the similarity to the two-mode approximation is lost at larger 
values of gi. 
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TABLE I. Intervals oi gi/oo for trapped oscillations 



Initial fl^Pl^^^ 




Q\/io (interval of trapped oscillations) 


0.1 


0.2 
0.4 


1.193 < gi/u < 2.522 
2.530 < gi/io < 4.635 


0.17 


0.2 
0.4 


1.982 < gi/oj < 3.154 
4.424 < gi/uj < 6.689 



IV. THE ATTRACTIVE NONLINEARITY 

In the version of the DWP model with the self-attraction, we start the simulations, 
without the application of the FRM, for values of the parameters similar to those of Ref.— : 
go = —1, gi = 0, Q = 0.1, Vq = 1, w = 0.5, and the initial value of (iV'P)max ~ 0.13. 
Results of the simulations are shown in Fig. |9l In this case, the unmanaged condensate is 
self-trapped in one well, oscillating in it at a frequency Wtrap.osc = 0.046. 
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FIG. 9. (Color online) The same as in Fig. [H but in the case of the attractive nonlinearity, with 
go = —1, gi = 0, and initial value (IV'P) . = 0.13. In the right panel, the solid and dashed 
lines represent the results obtained from the integration of the GPE and ODEs, respectively. This 
figure and others in this section correspond to the same parameters of the potential as above, i.e., 
n = 0.1, Vb = 1, and w = 0.5. 

In the case of the dominating attractive nonlinearity, the main effect of the FRM is 
opposite to that in the model with the repulsion: untrapping of the condensate which, in 
the unmanaged case, was trapped in the SSB state. For gi = 0.2, we have found that the 
untrapping occurs at 

0.1241 <uj < 0.1594 (14) 
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(or, almost exactly, Swtrap.osc < w < 4co'trap,osc)- This result is shown in Fig. [10] for uj = 0.13. 
In this case, the observed frequency of the oscillations of n{t) in the FRM- induced untrapped 
state is Wuntrap = 0.0318 ~ a;/4, suggesting that this may be a manifestation of a higher- 
order parametric resonance^. Notice that, as observed in the right panels of Figs. 191 and [TUl 
the agreement between the ODEs and the GPE results is good. 

We studied the oscillations of the untrapped state for other values of u from interval 
(fT4|) . It has been concluded that, at w = 0.12 and 0.14, the observed frequency of the 
untrapped oscillations again obeys the above-mentioned relation suggesting a higher-order 
parametric resonance: Wuntrap = 0.0298 ~ a;/4, and cuuntrap = 0.0346 ~ a;/4, respectively. 
In the left panel of Fig. (fTTj) . we present the relation between u and Wuntrap in the region 
of 0.124 < u < 0.159. For u< 0.124 and u >0.159, the oscillations of n{t) remain trapped 
in a single well, with the corresponding frequency almost the same as in the absence of the 
management, i.e., Wtrap.osc = 0.046. 
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FIG. 10. (Color online) The same as in Fig. [21 but in the model with the attractive nonlinearity, 
go = —1, with gi = 0.2, uj = 0.13, and initial value (I^P) = 0.13. In the right panel, the 
solid and dashed lines represent the results obtained from the integration of the GPE and ODEs, 
respectively. 

Next, we increased the value of the FRM amplitude, gi. For gi = 0.4, the untrapping 
occurred at 0.158 < u < 0.196 (which almost exactly means 4(X'trap,osc < ou < 5wtrap,osc)- As 
above, we examined the relation between u and the respective frequency of the untrapped 
oscillations, Wuntrap- The results verify the same conclusion as reached above, i.e., Wuntrap ~ 
u/i. In the right panel of Fig. [Til we present the values of u vs. cjuntrap in the region of 
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0.158 <uj < 0.196. 
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FIG. 11. (Color online) The relation between w and Wuntrap in the model with the attractive 
nonlinearity. Solid, and dotted lines correspond, respectively, to the result obtained from the 
GPE, and the estimation of Wuntrap = oj/A. Left panel: go = —1, gi = 0.2, 0.124 < a; < 0.159, and 
initial values = 0.13. Right panel: go = —1, gi = 0.4, 0.158 < cj < 0.196, and initial values 



As in the case of the repulsive nonUnearity, we have also examined the cases with gi > \go\, 
when the total nonlinearity coefficient periodically changes its sign. This was studied for 
Qi = 1.2 and gi = 1.5. In the former case, at 0.001 < u < 0.006 the population imbalance, 
n{t), oscillates in the untrapped state, see the left panels of Fig. [12] for u = 0.003. At 
0.007 < u < 0.014, we observe re-trapping after one or two untrapped oscillations, see the 
middle panels in Fig. [T2]for u = 0.007, while at a; > 0.014 the evolution of n{t) is irregular, 
as seen in the right panels of Fig. [T2]for u = 0.02. Similar results were obtained for gi = 1.5: 
in this case, n{t) oscillates without trapping at 0.001 < to < 0.004. For 0.004 < tu < 0.012, 
we observed re-trapping after one or two oscillations, and the evolution of n{t) is irregular 
at w > 0.012. As observed in the bottom panels of Fig. [TU it is obvious that the results 
obtained from the ODEs and the GPE are not in agreement for values of gi > 1, for similar 
reasons as considered above. 

V. CONCLUSIONS 

We have studied the trapping/detrapping of Bose-Einstein condensates confined in a 
DWP (double- well potential), under the action of the FRM (Feshbach-resonance-management) 
technique. Our model is based on the Gross-Pitaevskii equation with the time-dependent 



I^P = 0.35. 
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FIG. 12. (Color online) Spatiotemporal contour plots (top panels) and trajectories of n(i) (bot- 
tom panels) in the case of the attractive nonlinearity with qq = —1, gi = 1.2, and initial value 
(|^|2)^^^ = 0.13. Left panel: co = 0.003; middle panel: uj = 0.007; right panel: co = 0.02. Addi- 
tionally, In the bottom panels, the solid and dashed lines represent the results obtained from the 
integration of the GPE and ODEs, respectively. 

nonlinear coefficient whose average value corresponds to either repulsion or attraction. We 
have found that the two signs give rise to different effects, regarding the trapping/ detrapping 
of the condensate in the DWP. In the case of the self-repulsion, the BEC gets trapped in 
one of the two wells under the action of the FRM, if it was untrapped in the unmanaged 
state. We have found the range of values of the driving frequency leading to this effect. If 
the sign of the nonlinearity does not change due to the time-periodic modulation, and the 
densities are sufficiently low, as mentioned above, the frequency of the intrinsic oscillations 
of the FRM-induced trapped states is half the management frequency, which suggests a 
parametric resonance as an underlying mechanism. 

On the other hand, in the case of the self-attractive nonlinearity (on the average), we 
observe the opposite phenomenon: the application of the FRM results in untrapped oscil- 
lations, if the unmanaged condensate was self-trapped in one of the two wells. If the pe- 
riodically modulated nonlinearity does not change its sign, the frequency of the untrapped 
oscillations supported by the FRM is almost exactly equal to a quarter of the driving fre- 
quency, suggesting a higher-order parametric resonance as an underlying mechanism. It is 
worth mentioning here that, this frequency relation, is relevant, as in the previous case, for 
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sufficiently low densities. 

Our results were corroborated with ones derived in the framework of a semi-analytical 
approach. The latter, was based on an a Galerkin-type expansion that was used to describe 
the evolution of the wavefunctions at each well of the DWP. This led to a non-autonomous 
system of ODEs, which were solved numerically to find the relative difference in the atomic 
population between the two wells. In all cases, the ODE result was compared to the one ob- 
tained by the GPE model, and they were found to be in fairly good agreement: in particular, 
we found a quantitatively good agreement for small values of the modulation depth, gi < 1 
(in both repulsive and attractive cases), and a qualitatively good agreement for gi > I (for 
the repulsive case). Our approximation was found to fail only for (71 > 1 in which case the 
strong nonlinearity causes the involvement of additional modes in the dynamics. 

Detailed analysis of the resonances that may explain the character of the FRM-induced 
trapping and detrapping is a nontrivial problem, because the two-mode approximation, 
such as the one presented above or the one recently proposed in Ref.— , is insufficient for 
this purpose. The development of a more sophisticated finite-mode approximation, that 
may help uncover the resonances, is a subject for a separate study. Another challenging 
possibility is to consider similar dynamical effects in multidimensional settings. It may also 
be interesting to extend the analysis to the case of a binary BEG mixture trapped in a DWP. 
These problems will also be considered elsewhere. 
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